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COMBINATIONS  OF  GROUP  MEANS  SEWAl  RtCOBOS 

Abstract.  A  description  of  a  FORTRAN  IV  computer  pro- 
gram that  is  used  for  making  Scheffe's  S-contrasts  in  one-way 
analyses. 

Forestry  research  often  involves  experiments  in  which  many  treatments 
are  tested.  If  the  overall  test  of  treatment  equality  results  in  a  significant 
difference,  it  is  usually  necessary  to  make  further  tests  involving  linear 
combinations  of  the  treatment  responses.  When  many  a  posteriori  com- 
parisons are  necessary  and  treatments  have  unequal  sample  sizes,  the 
researcher  often  resorts  to  Scheffe's  S-contrasts  and  finds  himself  involved 
in  lengthy  calculations. 

Scheffe^  describes  the  method  of  making  contrasts  for  a  one-way  anal- 

I 

ysis.  In  short,  the  method  results  in  the  contrast  2  c.^(B.^  with  a  confidence 

1 

interval: 


I  c,P,  =1=  (M)(F«;  ,,^,)  (s)_ 


I 


Where: 


I 


Ci  =  a  contrast  multiplier  for  group  i,  and  2^=0. 

1 

=  arithmetic  mean  of  group  i. 


Scheffe,  Henry.  THE  ANALYSIS  OF  VARIANCE.  John  Wiley,  N.  Y.  477  pp.  1959. 


I    =  number  of  groups. 

F  =  Fisher's  F-value  for  I-l  and  N-I  degrees  of  freedom  at  the  chosen 

(a)  level  of  significance. 
N  =  total  sample  size  (all  groups  combined), 
s-  =  the  error  mean  square. 
J  J  =  sample  size  for  group  i. 

By  the  same  method,  we  may  calculate  an  F-value  for  an  individual 
contrast  as: 


(;■') 


(I.l)s^ 


and  we  may  check  the  F-value  against  a  tabulated  F  at  level  of  sig- 
nificance a. 

The  FORTRAN  IV  computer  program  described  here  is  simply  an 
automated  method  of  making  one-way  analysis  S-contrasts  for  linear 
combinations  of  group  means  with  constant  or  varying  sample  sizes. 

Description  of  Control  Deck 

The  user  must  supply  the  following  control  deck  as  data  to  be  oper- 
ated on  by  the  program: 

Card 

Card  No.  Columns  Content 

1  1-4       Number  of  treatments  or  groups, 

right  adjusted. 
5-8       Degrees  of  freedom  for  error, 
right  adjusted. 
9-17       Error  mean  square,  punched  with 

decimal  point. 
18-22       MEANS  punched   if  individual 
treatment  means  are  to  be  com- 
pared; blank  otherwise. 
2-1  1-10       Label  for  group. 

(1  for  each  group)        11-13       Sample  size,  right  adjusted 

14-20       Group  total,  punched  with  deci- 
mal point. 

21       Minus  sign  ( — )  if  multiplier  (c) 
for  contrast  is  negative. 
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22-25       Multiplier  (c)  for  contrast, 
punched  with  decimal  point. 

26-70  Repetition  of  format  for  columns 
21-25  for  each  additional  contrast 
to  be  made. 

The  program  is  limited  to: 

(1)  A  maximum  of  10  contrasts  (in  addition  to  the  contrasts  of  indi- 
vidual means). 

(2)  A  maximum  of  1,000  groups. 

Example 

If  we  assume  that  we  have  tested  the  effects  of  eight  fertilizers  on 
height  growth  of  red  pine  seedlings  and  that  the  treatments  are  repre- 
sented by  unequal  sample  sizes,  then  the  analysis  might  be: 


Source 

dj. 

SS 

MS. 

F. 

Fertilizer 

1 

94.0678 

13.438 

14.0 

Error 

134 

128.7060 

0.960 

Total 

141 

222.7738 

The  treatment  totals  and  their  sample  sizes  are: 

Total 

Treatment  Ji  ( groivth  in  height) 


1 

18 

36.7 

2 

19 

12.1 

3 

17 

11.8 

4 

15 

42.1 

5 

19 

38.9 

6 

20 

60.3 

7 

20 

40.2 

8 

14 

30.1 

8  134 

0. 

960MEANS 

TREAT. 

I 

18 

36.7 

.25  .25 

TREAT. 

2 

19 

12.1  .5 

-.5 

TREAT. 

3 

17 

11.8  .5 

-.5 

TREAT. 

4 

15 

42.1-.5 

-.5 

TREAT. 

5 

19 

38.9 

.25  .25 

TREAT. 

6 

20 

60. 3-. 5 

-.5 

TREAT. 

7 

20 

40.2 

.25  .25 

TREAT. 

8 

14 

30.1 

.25  .25 

Figure  1. — Control  deck. 


3 


S  CONTRASTS   INVOLVING  LINEAR  COMBINATIONS  OF  ARITHMETIC   MEANS  - 


EACH  CONTRAST  HAVING         7   AND     13i+   DEGREES   OF  FREEDOM 


TREAT . 

VS  . 

TREAT. 

2  1 

0.270AE 

01 

TREAT . 

VS. 

TREAT . 

3  1 

0.2353E 

01 

TREAT . 

VS  . 

TREAT . 

'i  1 

0.7177E 

00 

TREAT . 

VS  . 

TREAT . 

0.9890E-04 

TREAT . 

VS  . 

TREAT . 

6  1 

0.13A3E 

01 

TREAT. 

VS  . 

TREAT  . 

7  1 

0.1177E- 

-02 

TREAT  , 

VS. 

TREAT. 

8  1 

0.1A47E-01 

TREAT . 

2 

VS  . 

TREAT. 

3  1 

0.4380E- 

-02 

TREAT . 

2 

vs. 

TREAT . 

A  1 

0.5873E 

01 

TREAT . 

2 

VS . 

TREAT. 

5 

0.2813E 

01 

TREAT • 

2 

vs. 

TR  EAT . 

6  1 

0.8200E 

01 

TREAT . 

2 

vs. 

TREAT. 

7  1 

0.2734E 

01 

TREAT  . 

2 

VS  . 

TREAT . 

8 

0.2746E 

01 

TREAT . 

3 

VS . 

TR  E AT . 

A 

0.5292E 

01 

TREAT. 

3 

VS. 

TREAT. 

5 

0 . 244  5E 

0  1 

TREAT. 

3 

vs. 

TREAT . 

6 

0.7366E 

01 

TREAT . 

3 

VS. 

TREAT . 

7 

0.2368E 

01 

TREAT . 

3 

vs. 

TREAT. 

8 

0.2422E 

01 

TREAT . 

4 

vs. 

TREAT . 

5 

0.7I92E 

00 

TREAT . 

vs. 

TREAT, 

6 

0.5536E- 

-01 

TREAT . 

<t 

VS  . 

TREAT . 

7 

0.8095E 

00 

TREAT. 

vs. 

TREAT. 

e 

0.4647E 

00 

TREAT. 

5 

VS. 

TREA1  . 

6 

0.1358E 

01 

TREAT. 

5 

VS. 

TREAT. 

7 

0.2025E- 

-02 

TREAT. 

5 

VS. 

TREAT. 

8 

0.  1263E- 

-01 

TREAT. 

6 

vs. 

TREAT. 

7 

0.1503E 

01 

TREAT. 

6 

VS. 

TREAT. 

8 

0.9169E 

00 

TREAT. 

7 

vs. 

TREAT. 

8 

0.2402E- 

-01 

0.5000E 

00 

X 

TREAT. 

2 

0-5000E 

00 

X 

TREAT. 

3 

-0.5000E 

00 

X 

TREAT. 

4 

-0.5000E 

00 

X 

TREAT. 

6 
F 

0.25GOE 

00 

X 

TREAT. 

1 

-0.5000E 

00 

X 

TREAT. 

2 

-0.5000E 

00 

X 

TREAT. 

3 

0.2500E 

00 

X 

TREAT. 

5 

0.2500E 

00 

X 

TREAT. 

7 

0.2500E 

00 

X 

TREAT. 

8 
F 

0.2500E 

00 

X 

TREAT. 

1 

-0.5000E 

00 

X 

TREAT. 

4 

0.2500E 

00 

X 

TREAT. 

5 

-0.5000E 

00 

X 

TREAT. 

6 

0.2500E 

00 

X 

TREAT. 

7 

0.2500E 

00 

X 

TREAT. 

8 

F=  0.2466E  01 

Figure  2. — Program  output 
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SIBFTC  STEST 

DIMENSION  LABEL ( I 000 , 3 ) , X J ( 1 000 ) , TOT ( 1000 ) , C (I 000 , 10 ) 
READ<  5.1) IGRPS, I  OF , SQUARE , MEAN, S 

1  FORMAT(2IA,F9.0»A<»,A1) 
DATA  MtA,SN/4HMEAN, IHS/ 
XGRPS=IGRPS 

CF=IDF 

IGRPDF=IGRPS-1 
DO  2   1=1, IGRPS 

2  READ( 5,3) (LABEL ( I , J ) , J= 1 , 3 ) , X J ( I ) ,TOT( I ) , ( C ( I , J ) , J= 1 , 10 ) 

3  FORMAT(2A4,A2,F3.0.F7.0, 10F5.0) 
kRITEC 6,4) IGRPDF, IDF  • 

4  F0RMAT(85H1S  CONTRASTS   INVOLVING   LINEAR  COMBINATIONS  OF  ARITHMETIC 
1   MEANS  -  EACH  CONTRAST  HAVINGI5,4H  ANDI5,20H  DEGREES  OF   FREEDOM  ) 

IF(MEA.NE.MEAN.OR.SN.NE.S)GO  TO  7 

CO  5   1=1, IGRPDF 

IZ=I*1 

CO  5  K=IZ, IGRPS 

CIFF=TOT( I )/XJ( I )-TOT(K )/XJtK) 
SUM=1./XJ( I )+l./XJ(K) 

F=(DIFF**2)/( (XGRPS-1. )*SQUARE*SUM) 

5  WRITE(6,6)(LABEL( I,J),J=1,3), (LABEL(K,J),J=1,3) ,F 

6  FnRMAT(2H0   2A4,A2,5H  VS.    2A4,A2,4H,  F=E11.4) 

7  CO  13  J=l,10 
V«RITE(6.8) 

8  FORMATdH  ) 

CO  9  1=1, IGRPS 
IF(C(I,J).NE.0.0)G0  TO  10 

9  CONTINUE 
GO  TO  15 

10  CIFF=0. 
SUM=0. 

CO  12   1=1, IGRPS 

CIFF=DIFF+C( I,J)*(T0T(1 )/XJ(I ) ) 

IF(C{I,J).NE.0.)WRITE(6,11)CII,J) , ( L ABEL ( I , K ) , K= 1 , 3  ) 

11  FORMATdH  E11.4,3H  X  2A4,A2) 

12  SUM=SUM+IC( I,J)**2)/XJ( I ) 
F={DIFF**2)/( (XGRPS-1. ) *SQUARE*SUM ) 

13  ViRITE(6,  14)F 

14  FORMATdH  2 IX  ,  2HF  =  E  1 1 .  4 ) 

15  STOP 
END 


Figure  3. — Program  listing. 


Now,  suppose  that  we  want  to  make  contrasts  of: 

(1)  individual  treatment  response  means. 

(2)  average  of  responses  to  treatments  2  and  3  versus  average  of  4 
and  5. 

(3)  average  of  responses  to  treatments  1,  3,  7,  and  8  versus  average 
of  2  and  3. 

(4)  average  of  responses  to  treatments  1,  5,  7,  and  8  versus  average 
of  4  and  6. 
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Then  the  control  deck  for  this  job  would  be  that  shown  in  figure  1.  The 
results  are  shown  in  figure  2.  If  we  check  the  calculated  F-values^ 
against  tabulated  Fj  ^gj^  =  2.08,  we  find  the  following  contrasts  to 
be  significant: 

(1)  Treatments  2  and  3,  taken  separately,  versus  each  of  the  remaining 
treatment  means  (1,4-8). 

(2)  Contrasts  (2),  (3),  and  (4),  as  stated  above. 

Figure  3  contains  a  program  listing. 

—  WARREN  E.  FRAYER3 

Assistant  Professor  of  Forest  Mensuration 
Colorado  State  University 
Fort  Collins,  Colorado 


'  Note  that  the  floating-point  multipliers  and  F- values  contain  a  decimal  value,  and,  after 
the  E,  the  number  of  places  to  the  left  ( — )  or  right  that  the  decimal  point  must  be  moved. 


'  When  this  paper  was  prepared,  Dr.  Frayer  was  a  research  forester  with  the  Northeastern 
Forest  Experiment  Station,  Forest  Service,  U.  S.  Department  of  Agriculture,  Upper  Darby,  Pa. 


